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Materials which can undergo slow diffusive transformations as well as fast displacive transfor- 
mations are studied using the phase- field method. The model captures the essential features of 
the time-temperature-transformation (TTT) diagrams, continuous cooling transformation (CCT) 
diagrams, and microstructure formation of these alloys. In some materials systems there can exist an 
intrinsic volume change associated with these transformations. We show that these coherency strains 
can stabilize mixed microstructures (such as retained austenite-martensite and pearlite-martensite 
mixtures) by an interplay between diffusive and displacive mechanisms, which can alter TTT and 
CCT diagrams. Depending on the conditions there can be competitive or cooperative nucleation of 
the two kinds of phases. The model also shows that small differences in volume changes can have 
noticeable effects on the early stages of martensite formation and on the resulting microstructures. 
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I. INTRODUCTION 

Solid-to-solid phase transformations have been tra- 
ditionally classified as diffusive and displacive, de- 
pending on the underlying kinetics.^ Diffusive phase 
transformations, such as spinodal decomposition and 
precipitation in alloys, are slow because they require 
long-range motion of the atoms. ^ In displacive phase 
transformations, on the other hand, the crystal structure 
changes through a unit cell distortion. Since atoms 
move over very short distances, these transformations 
can be very rapid. First order displacive transformations 
— also referred to as martensitic transformations — are 
responsible for the shape-memory effect and for the pseu- 
doelastic behavior ^^i^ They are usually accompanied by 
a spontaneous strain. The high temperature cubic phase 
(austenite) transforms to a low temperature tetragonal, 
orthorhombic or monoclinic phase (martensite).^ 

Diffusive and displacive phase transformations may 
interact or compete with each other. The best known 
material system where this interplay is observed is 
eutectoid steel. The high temperature austenite phase 
typically decomposes into pearlite, i.e. ferrite plus 
cementite (iron carbide, FeaC), by a diffusive process. 
For fast cooling rates, a transformation from austenite 
(face-centered cubic) to metastable martensite (base- 
centered tetragonal) may instead take place. 

Interplay between diffusive and displacive transforma- 
tions has also been observed in other materials systems, 
for instance Ti-Al-Nb^ Cu-Al-Ag^ Cu-Zn-Al^^ and 
Pu-Ga alloysi^ In NiTi shape-memory alloys, the forma- 
tion of precipitates by a diffusive phase transformation 
influences the mechanical response and the stress-induced 
martensitic transformations!^ 

Information on diffusive and displacive phase trans- 
formations and their kinetics can be conveniently repre- 



sented on time-transformation-temperature (TTT) and 
continuous cooling transformation (CCT) diagrams: the 
volume fractions of the different phases or microstruc- 
tures are plotted as a function of time at different 
temperatures. These diagrams account for metastable 
phases such as martensite which do not appear in phase 
diagrams. 

In order to describe the full microstructural complexity 
generated by the transformations, it is necessary to 
use numerical, rather than analytical, treatments. The 
phase-field method (also known as the time- dependent 
Ginzburg-Landau method) is a powerful technique to de- 
scribe the microstructures and kinetics of phase transfor- 
mations. In particular, one needs not make any assump- 
tion regarding the arrangement of the different phases: 
microstructures are an output of the simulations. It has 
been extensively used to study microstructural evolution 
in diffusiveiaiiO^OS and martensitioii^^^^^^^^ 
transformations. However, apart from one work on 
a TTT diagram for ferrite and martensite in iron^ 
diffusive and displacive transformations have never been 
studied together using the phase-field method. 

In a previous work^^ we introduced a phase-field model 
to study the interplay between diffusive and displacive 
transformations. TTT diagrams and microstructures 
were obtained. If the atomic volume of martensite is 
different from that of austenite or if there is a lattice 
mismatch due to phase separation, phase transformations 
create hydrostatic strain. An important conclusion 
from our previous work was that these intrinsic volume 
changes can stabilize mixed microstructures. In the 
present article we systematically study the effect of such 
strains on microstructures, TTT diagrams, and CCT 
diagrams. 

Section [III presents a phase-field model suitable for 
a system which can undergo phase separation as well 
as a square-to-rectangle martensitic transformation. In 
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Sec. Uni the model is used to obtain TTT diagrams 
and microstructures in the absence of volume change. 
Sections |lVl and [Vl focus on the effect on isothermal trans- 
formations of a volume change associated with martensite 
and pearlite formations, respectively. Section lVll presents 
continuous cooling results. 



II. PHASE-FIELD MODEL 

Since the system can undergo both a phase separation 
and a martensitic transformation, we consider three 
'phases': austenite, martensite (of which there exists 
two variants, elongated along [01] and along [10], 
respectively), and pearlite (of which there are two 
kinds, with compositions higher and lower than that 
of austenite). We therefore label 'pearlite' any phase- 
separated region, even if the microstructure does not 
correspond to what one would typically call pearlite. 
Further, although the vocabulary used to designate 
phases may be reminiscent of steel, this is simply to avoid 
cumbersome periphrases. The model presented here is 
not specific to steel, neither is it claimed that it is a 
realistic representation of steel. 



A. Free energy 

We use two kinds of variables: composition and strains, 
c is the composition, it is conserved. Since only its 
variations are relevant to diffusion, c can be defined to a 
constant: we set c = as the composition of austenite 
and pearlite corresponds to c = ±co. ei, 62, and 63 are 
non-conserved variables, ei is the hydrostatic strain, 62 
is the deviatoric strain, and 63 is the shear strain: 
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The {sij} are the linearized strain tensor components. 
The free energy of the system is expressed as^^ 

G = j (sei + 5ch + 5cpi + y II Vcf + ^ II Veaf^dr, 

where g^i is the non-linear elastic free energy density for 
a square-to-rectangle martensitic transition : ^^i-*^^ 
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Here T is the dimensionless temperature and Tm and Tp 
are constants pertaining to the austenite-martensite and 
austenite-pearlite phase transformations respectively. 
The parameters X12 and xic are related to the volumetric 
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FIG. 1: The bulk energy as a function of 62 and c at (a) T — 
0.2, (b) T = 0.7, and (c) T = 1.5. Darker areas correspond 
to lower energies. 



strains coming from martensite and pearlite formations, 
respectively.^^ 

The chemical free energy is given byi^ 



B2 T— Tp 2 , ^4 4 



(3) 



and a coupling between elastic distortions and composi- 
tion is introduced as 



^cpl = X2c 6^(62)^ 



(4) 



Since we are interested in a qualitative understanding 
of the physical mechanisms we choose simple values for 
the parameters. Yet we ensure that the energy of pearlite 
is always lower than that of martensite and that pearlite 
formation is slower than martensitic transformation by 
orders of magnitude. Ai = 1, ^22 = 2, A2A = 4, 
A2Q = 9.6, A3 = 1, ^2 = 6, ^4 = 12, X2c = 5, Tm = 0.5, 
Tp = 1, = 2, and ke2 = 0.1. The coupling constants 
X12 and xic are varied in different cases to understand 
the effect of volume changes: in Sec. IIIII both are set to 
zero. Sec. IIVI focuses on X12 7^ and Sec.[Vlon xic 7^ 0. 

The homogeneous part of the free energy is depicted 
in Fig. [1] as a function of 62 and c at different 
temperatures. The different phases can be identified as 
follows: austenite corresponds to c = and 62 = 0, 
martensite to c = and 62 7^ 0, and pearlite corresponds 
to c 7^ and 62 = 0. Above Tp only austenite is stable. 
Fig. l(c)[ Between Tm and Tp austenite and martensite 
are unstable and pearlite is the ground state. Fig. |l(b)[ 
Below Tm , pearl ite is the ground state and martensite is 
metastable. Fig. |l(a)[ 



B. Evolution equations 

The evolution of the composition is described by the 
Cahn-Hillard equationi^ 



dc{r,t) . SG 



dt 



(5) 



where 5 is the functional derivative and M is the 
temperature-dependent mobihty, 



M = Mo exp(-g/T). 



(6) 
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FIG. 2: Time-temperature-transformation diagrams for several values of xi2 and xic. A: at least 10% austenite, M: at least 
10% martensite, and P: at least 10% pearlite. Points a-m correspond to the conditions at which the microstructures in Figs.[3l 
m [11] and [12] are obtained. 



This temperature dependence of the mobility critically 
influences the TTT diagrams. 

The evolution of the displacements is described by2^ 



-^'"^('■'')=^^^|M+,vV(r,t), (7a) 
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5 Eij (r, t) ' 
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where p is a density, {dij} are stresses, and v is the 
time derivative of the displacements, u. The second term 
on the right-hand side in Eq. ([Taj) is a viscous damping 
term; it is a simpliflcation of the more general damping 
of Ref. m. 

After the displacements are calculated from Eq. ([7j), 
the strains are obtained as their derivatives. Then ei, 62, 
and 63 are obtained via Eq. ^ and are used in Eqs. ([2]), 
(|3|), and (|4]) to calculate the energy. The energy will in 
turn be used to evolve the displacements through Eq. (O 
at the next time step. 

All simulations are two-dimensional (128 x 128 lattice) 
with periodic boundaries A flnite difference scheme 
with Ax = 1 and At = 0.2 is used. We use r] = 0.01, 



p = 1, Mo = 2, and Q = 5. The initial system, made 
of 100% austenite, includes random fluctuations around 
c = and u = 0. The system is quenched instantaneously 
to temperature T and held at this temperature. For each 
value of T we record the times at which 10% martensite 
and 10% pearlite form, as well as the time at which the 
austenite content drops below 10%. 



III. IN THE ABSENCE OF VOLUME CHANGE 

In this section we consider the case of X12 = and 
xic = 0, i.e. that there is no hydrostatic strain associated 
with either martensite or pearlite formation. 



A. TTT diagram 



Figure |2(a)| shows the resulting TTT diagram. The 
austenite-pearlite phase transformation requires diffu- 
sion and therefore time. At low temperature, diffusion 
is slow and so is pearlite formation. At temperatures 
close to Tp = 1 the driving force is small and pearlite 




FIG. 3: (color online) Pear lite microstructures (light blue: 
'cementite' and dark blue: 'ferrite'). (a): xi2 = xic = at 
T = 0.5 and t = 2 000, point a in Fig. |2(a)| (b): 
at T = 0.9 and t = 2 000, point b in Fig. |2(a) 
xic = 1 at T = 0.85 and t = 2 000, point f in Fig. [2(c)] (d): 
xi2 = xic = 1 at T = 0.8 and t = 10 000, point g in Fig. ]2(d)| 
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FIG. 4: (color online) Microstructures for xi2 = xic = at 
T = 0.49 and (a) t = 50, (b) t = 4 000, and (c) t = 12 000. 
They correspond to points c, d, and e in Fig. 2(a) Red and 
yellow: martensite; light and dark blue: pearlite. 



formation again is slow. Consequently there exists an 
intermediate temperature at which pearlite formation is 
the fastest, this accounts for the C-curve around T = 0.8 
m Fig. [2(a)| This is a typical feature of experimental 
TTT diagrams. 

Above Tm = 0.5, only pearlite is stable. Therefore it 
will necessarily form, albeit slowly. Below Tm, although 
pearlite is still the ground state, martensite is metastable. 
Unlike pearlite, martensite forms through a displacive 
mechanism, which does not require long-range motion of 
atoms. The phase transformations are then controlled 
by kinetics and the fast martensitic transformation takes 
place instead of the thermodynamically favorable (but 
slow) pearlite formation: Fig. 2(a) shows that below 
T ^ 0.5 (the 'martensite start temperature') austenite 
transforms to martensite. Notice that eventually pearlite 
will form even below T = 0.5 as it is the ground 
state. Pearlite formation below the martensite start 
temperature has been observed experimentally*2i 



B. Microstructure evolution 



2(a)| shows the existence 



The TTT diagram in Fig. 
of various microstructures at different temperatures and 
times. The typical length scale of pearlite depends on 
the diffusion length and therefore on temperature: at 
low temperatures pearlite is fi ne — Fig. 3(a) — and it 



is coarser close to Tp — Fig. |3(b)| — , consistent with 
experimental observations. Note that the amounts of the 
two components of pearlite are equal because, as Fig. [1] 
shows, their compositions are symmetric with respect to 
that of austenite. 

Figure 2] shows the microstructure evolution at T = 
0.49 (just below the martensite start temperature). 
At t = 50 only martensite is found. Fig. |4(a)[ At 
t = 4 000, pearlite has already started forming at the 
interfaces between the martensite variants aligned along 
[1 1] and those along [1 1], Fig. |4(b)] that is at interfaces 
which are relatively high in energy. For longer times 
pearlite grows at the expense of martensite, which finally 
disappears [Fig. 4(c)| . This appearance of pearlite in the 
martensite region occurs at progressively later times at 
lower temperatures. Figure 4(c) (obtained at T = 0.49) 



is similar to Fig. 3(a), which is obtained at T = 0.5. 
However, the two microstructures were the outcomes of 
very different routes. 

As discussed earlier, below T = 0.5, pearlite is the 
ground state and martensite is only metastable. Thus, 
the martensite forms only because pearlite formation 
is slow at this temperature. It is therefore expected 
that pearlite can eventually form in a martensitic system 
as this will decrease the energy. However there exists 
an energy barrier because — unlike austenite which is 
unstable below T = 1 — martensite is metastable. This 
makes pearlite nucleation off martensite slower than 
pearlite nucleation off austenite and accounts for the 
discontinuity of the 10% pearlite hue across T = 0.5 in 
Fig. [2(^01^-1 



IV. EFFECT OF MARTENSITE VOLUME 
CHANGE (xi2 / 0) 

In the previous section, we assumed that there was no 
strain associated with pearlite or martensite formation, 
i.e. Xi2 = Xlc = in Eq. ([2]). We now relax this 
constraint and look at the effect of the term xi2 (^2)^ in 
the energy expression, keeping xic set to zero. This term 
couples deviatoric strain (i.e. martensite) and hydrostatic 
strain: X12 7^ associates a net volume change with the 
martensitic transformation. 



A. TTT diagram 



Figure [2(b)] shows the TTT dia gram for X12 = 1 and 
Xlc = 0. There are two noticeable differences compared 
to the TTT diagram obtained with X12 = xic = 0, 
Fig. 2(a): there exists a martensite start temperature 



as well as a martensite finish temperature and there is 
hardly any discontinuity of the 10% pearlite line across 
T = 0.5. 

The martensite start and finish temperatures are due 
to the presence of retained austenite between T ^ 0.45 
and T = 0.5. The hydrostatic strain associated with the 
martensitic transformation results in coherency stresses. 
Therefore martensite growth in the austenite matrix is 
arrested before it is complete. The martensite/austenite 
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FIG. 5: (color online) Microstructures at T = 0.49 and t = 50. 
From left to right: xi2 = 0, xi2 = 0.85, xi2 = 1, and 
xi2 = 1.5. Top: xic = and bottom: xic = 1. Red and 
yellow: martensite; light and dark blue: pearlite. 



ratio increases for decreasing temperature and the system 
is purely martensitic at very low temperature. This 
trend is consistent with experimental results: in Cu- 
17.0% Zn-13.7 at.% Al, the volume change is small 
and the difference between martensite start and finish 
temperatures is 30 °C,^^ whereas it is around 100 °C in 
steel, where the volume change is larger. 

When xi2 = xic = 0, below T = 0.5 pearlite 
can nucleate heterogeneously at the interface between 
martensite variants (Fig. 2]); this process is rather slow. 
When Xi2 = 1, even below T = 0.5 some austenite 
remains and pearlite formation proceeds in this retained 
austenite. Since austenite is unstable there is no 
nucleation barrier for pearlite formation and there is 
hardly any discontinuity of the 10% pearlite line. This 
is unlike in Fig. |4(b)| where pearlite must nucleate at 
the boundary between martensite variants as there is no 
retained austenite. The pearlite nucleation is slower at 
lower T only because the diffusivity is lower. 



B. Microstructures 

The effect of Xi2 on the microstructure at T = 0.49 

If Xi2 7^ 0, there is 



can be seen in Figs. 5(a)-5(d) 



a competition between the driving force favoring the 
transformation and the coherency stress arising out of 
the volume change: the larger Xi2, the larger the stress 
and the lower the equilibrium martensite fraction. This 
gives rise to several microstructures for increasing xi2. If 
Xi2 is large the system is mostly austenitic, Fig. 5(d) At 
intermediate values the system looks like a checkerboard 
of alternating martensite and austenite, Fig. 5(c)| (this 
kind of microstructure has been observed experimentally 
by Le Bouar et al.^^). At lower xi2 the system is mostly 
martensite with small austenitic regions. Fig. |5(b)[ For 
even lower values of Xi2 the stress is not sufficient to 
make a somewhat relaxed austenite-martensite mixture 
more favorable than stressed martensite and the system 
is made of pure martensite. 




FIG. 6: (color online) Microstructures for xi2 = 0.85 and 
xic = at T = 0.49: (a) t = 1 000, (h) t = 2 000, and (c) 
t = 3 500. Red and yellow: martensite; light and dark blue: 
pearlite. 



Comparison of Figs. [S^e} - 5(h) to Figs. |5 (a)"] - |5 (d) | shows 
that, at T = 0.49 and t = 50, the microstructures do not 
depend on the value of xic nearly as much as on the value 
of X12: the mixed austenite-martensite microstructures 
are controlled by the value of X12. This is because xic is 
important only in the presence of pearlite. 

Along with Fig. |5(b)[ Fig. [6] shows the microstructure 
evolution at T = 0.49 for X12 = 0.85 and xic = 0. Unlike 
pearlite nucleated from austenite, Fig. |3(a)| pearlite here 
exhibits some texture along [1 1], Fig. |6(c)[ This texture 
is due to the way pearlite grew, not to thermodynamics; 
it is therefore different from Figs. 3(c) and |3(d)| (see 
This alignment of pearl- 



ite perpendicularly 



Sec. IVBl) 

to the interface where it nucleated has been extensively 
described. 



C. Interface orientation 

As expected from elastic compatibility, martensite- 
martensite interfaces are oriented along (11), as seen 
However, the orientation of martensite- 



in Fig. 5(a) 



austenite interfaces is different, as shown in Figs. 5(b) 
|5(d)| (the same holds for martensite-pearlite interfaces). 

We consider an interface (martensite-martensite, 
martensite-pearlite, or martensite-austenite) at an angle 
with the y-axis and determine what values of are 
allowed by elastic compatibility, ei, 62, and 63 depend 
only on x cos 6 -\-ysm6 and, for i = 1, 2, 3, 



dx'^ 



e'' cos and 



e'' sin^ ( 



(8) 



As the only contribution of 63 to the free energy is 
through the (63)^ term, energy minimization with respect 
to 63 gives 63 = 0. Elastic compatibility then requires^ 



d'^ei d^e\ d^e2 d^e2 



0. 



(9) 



dx^ dy^ dx^ dy^ 

If x\cC + x\2 (62)^ = in Eq. ([2]), energy minimization 
gives ei = 0. This occurs if (i) x\2 = for a 
martensite-martensite or a martensite-austenite inter- 
face, (ii) x\c = for a pearlite-austenite interface, or (iii) 
^12 = ^ic = 0- Equations (O and ([9]) then give 



cos26> = 0. 
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FIG. 7: (color online) Microstructures at T = 0.49 and 
xic = for various times and values of xi2. Red and yellow: 
martensite; green: austenite. 



In order to satisfy compatibility, interfaces must 
be aligned along (11), as observed for martensite- 
martensite interfaces in Figs. 5(a) and |5(e) for instance. 

If xicC + Xi2 (62)^ is not equal to zero, neither is the 
equilibrium value of ei. Each of the two phases may 
then have a different value of ei and 62- Calling Sei and 
6e2 the jumps of ei and 62 across the interface, Eqs. ([8]) 
and (O give 



Sei cos + 6ei sin — 6e2 cos + 6e2 sin = 0. 



This condition is satisfied if 



062 



(11) 



For an interface between two martensite variants, ei is 
constant {Sei = 0) and one recovers Eq. (p!Q|) . For 
martensite-austenite and martensite-pearlite interfaces 
dei needs not be zero. In the case of Fig. 5(g), for 
instance, the values of ei and 62 give 6 ^ ±30° or 6 ^ 
±60°, which is consistent with the interface orientations 
observed in the figure. 

In the case of pearlite-pearlite interfaces (i.e. ferrite- 
cementite interfaces), 62 = everywhere and only ei 
varies spatially if Xic 7^ 0. Equation (pT]) then gives 



FIG. 8: (color online) Microstructures at T = 0.49 and xic = 
for X12 = 0.75 (top) and 0:12 = 0.85 (bottom). Red and 
yellow: martensite; green: austenite. 



cos26> — 0, i.e. the interfaces are along (11), as seen 



3(d)[ Although this result is identical 



in Figs. |3(c)| and 

to Eq. (flop / the latter does not apply to pearlite as it 
assumes that ei = and 62 7^ 0, which is the opposite of 
the pearlite case. If xic = then both ei and 62 are zero 
everywhere and there is no constraint on the orientation 
of the interfaces, as shown by Figs. |3(a")] and 3(b) 



D. Early stages of martensite formation 



Figures 5(a) - 5(d) show four different microstructures: 
pure martensite, mostly martensitic with small austenitic 
regions, patterned (checkerboard) martensite-austenite 
mixture, and isolated martensite grains in an austenite 
matrix. Figure [71 shows the early stages of martensite 
formation leading to these microstructures.^^ 

The initial stages of martensite formation are very 
similar for all values of X12 and the differentiation 
between these various types of microstructures occurs 
rather late. At t = 0.85 and t = 0.9 the differences 
between X12 = 0.75, X12 = 0.85, and X12 = 1 are very 
minor. Differentiation occurs at t = 1: continuous strips 
form at X12 = 0.75 and ^12 = 0.85 [Figs. |7(g)| and |7(k)| 
but not at X12 = 1, Fig. |7(o)[ 

Differentiation between X12 = 0.75 and X12 = 0.85 
occurs even later (see Fig. [8]). Not only will the final 
martensite fraction be different but the orientation too 
will differ. Att = 1.5, X12 = 0.75 and X12 = 0.85 are very 
similar. Figs. [7(h) and |7(l)[ but for X12 = 0.75 martensite 
will finally align along [1 IJ , Fig. |8(d)] whereas in the case 
of X12 =0.85 some auste nite is retained and martensite 
is aligned along [1 1], Fig. |8(h)[ 

If X12 is smaller than about 0.8, the system will 
transform to pure martensite. Up to X12 ~ 0.65, the 
martensitic transformation completes very rapidly after 
it starts (by t = 1). For larger values of X12 (e.g. X12 = 
0.75) there is a two-step process: first 80-90% martensite 
forms, as in Fig. |7(h)[ and the remaining austenite 
transforms to martensite only much later. Fig. |8(d)[ 




At t = 0.85, for 0.75 < Xi2 < 1.5 small martensite 
grains are aligned along (11). When xi2 = 1-5 this 
alignment is lost at later times (unlike what can be 
observed at lower values of X12) but the shape of 
the grains is roughly conserved [compare Fig. |5(d)| to 
Fig.[7(t)]. 

In spite of these differences, several features seem 
independent of the value of X12: (i) the transformation 
starts around t = 0.85, (ii) the final martensite fraction 
is reached very quickly (except for X12 = 0.75), and 
(iii) there is some form of alignment along (11) at the 
initial stages. Although all martensite formation occurs 
between t = 0.85 and t = 1, there may be non- negligible 
microstructural evolution up to t ^ 12, i.e. there can be 
two steps: first martensite forms very fast (within about 
0.15 time units) and then its structure evolves at constant 
martensite fraction (which can take up to 10 time units). 



V. EFFECT OF PEARLITE LATTICE 
MISMATCH (xie / 0) 

The term xic c in Eq. ([2]) corresponds to a lattice 
mismatch between ferrite and cementite: one expands 
and the other contracts; these strains cancel out and 
there is no net volume change associated with the 
austenite-pearlite transformation. Nevertheless, the co- 
herency strains generated by this volume change crucially 
influence the microstructures and TTT diagrams. 



A. TTT diagram 



Figure [2(c)] shows the TTT dia gram corresponding to 
X12 = and xic = 1. The main difference from the 
TTT diagram obtained with X12 = xic = 0, Fig. 2(a)[ is 
that martensite can be found above T = 0.5. There is 
then a possibility to observe a coexistence of martensite 



and pearlite. Indeed, due to the coherency strains and 
elastic compatibility, a mixture of pearlite and martensite 
is more stable than pearlite on its own (this remains true 
up to T ^ 0.7). 

This splitting between pure pearlite and pearlite- 
martensite mixture of what is the pearlite region in 
Fig. 2(a) is akin to the pearlite-bainite transition in 



steel. However, the corresponding microstructure, shown 
in Fig. |9(e) , is clearly different from that of bainite: 
in bainite structural distortion and diffusion occur in 
the same domain, whereas in the microstructure shown 
here diffusive and displacive mechanisms are involved in 
physically distinct regions. 



B. Microstructures 

Figure [9] shows the effect of xic on the microstructure. 
As mentioned in Se c. HVl X12 helps to retain austenite 
at short times: Figs. 5(a) and 5(e)| (which corresp ond to 
X12 = 0) show pure martensite whereas in Figs. |5(b)| - 
|5(d)| and 5(f)| - |5(h)] (corresponding to non-zero X12) there 
is a mixture of martensite and austenite. The long-time 
microstructures on the other hand are controlled by the 
value of xic'. in Figs. |9(a)| and 9(g)| xi c = and there 
is only pearlite whereas Figs. 9(b)| - [9(f)" and 9(h)f 9(1)] — 
which correspond to Xic > — show both pearlite and 
martensite. These pearlite-martensite microstructures 
do not disappear in longer simulations, xic stabilizes 
martensite at long times for a similar reason why X12 
stabilizes austenite at short times: when xic = 1 a purely 
pearlitic system is unstable as the lattice mismatch 
strains generated by the phase separation tend to cause 
a stress-induced martensitic transformation. 

Mixed microstructures can exist up to T ^ 0.7. 
Figure [To] shows that at T = 0.68, the pearlite-martensite 
mixtures are more regular than at lower temperature. 
Fig. |9(e) and |9(k)[ The greater diffusion allows the 
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FIG. 10: (color online) Microstructure at T = 0.68 and xic = 
1. (a): xi2 = and t = 40 000, (b): xi2 = 1 and t = 200 000. 




FIG. 11: (color online) Microstructures for xi2 = xic = 1 at 
T = 0.49: (a) t = 1 000, (b) t = 1 500, and (c) t = 2 000. 
They correspond to points k, 1, and m in Fig. 2(d) Red and 
yellow: martensite; light and dark blue: pearlite. 



systerQ to evolve towards lower Daininia, i.e. raore regular 
states. However, one should also note that this kind of 
self-assembly is due in part to the periodic boundaries 
which allow only discrete periods (128, 64, etc.). 

Figures 3(c)| and |3(d)| show that if xic = 1, pearlite 
is textured [compare to Fig. |3(b)| . If xic 7^ 0, one 
component of pearlite has a positive value of ei and 
the other a negative value. Elastic compatibility then 
dictates that the interfaces be along (1 1) as it is the case 
for martensite (see Sec. II V C]) . 



C. Both X12 and xic / 



Figure [2(d)] shows the TTT dia gram for X12 = xic = 1. 
It includes the features from the TTT diagrams obtained 
with X12 = 1 and xic = [Fig. |2(b)| , and with X12 = 
and xic = 1 [Fig. |2(c)| : there exist martensite start 
and finish temperatures, martensite can be found above 
T = 0.5, and there is no discontinuity across T = 0.5. 

Figures 5(g), |ll(a)| - p~l(c) , and |9(k)| show the time 
evolution of the microstructure at T = 0.49 for X12 = 
xic = 1. By t = 1 000, pearlite nucleated in the 
retained austenite and new martensite formed where 
there was none before (pure pearlite cannot form due 



to the lattice mismatch stresses it generates). Fig. 11(a) 



After austenite disappears pearlite grows at the expense 
of the large grains of 'primary' martensite. Fig. 11(b) 



The microstructure 



which finally disappear. Fig. |ll(c) 
then coarsens. Fig. |9(k)[ 

From large martensite grains alternating with areas 
completely devoid of martensite at t = 50 [Fig. |5(g)| , 
the system evolves to a state where martensite is more 
homogeneously distributed [Fig. ll(c)| through a double 




FIG. 12: (color online) Microstructures at T = 0.6 for X12 = 
xic = 1. (a) t = 50, (b): t = 100, a nd (c) : t = 200. They 
correspond to points h, i, and j in Fig. 2(d) Red and yellow: 
martensite; light and dark blue: pearlite; green: austenite. 



mechanism of martensite formation and martensite 
destruction, and finally to a coarsened microstructure. 
Fig. |9(k)[ The feature size thus goes from large to small 
to medium with increasing time. Consequently the 'final' 
structure is independent of the initial one, both in terms 
of grain size and interface orientation [compare Fig. |9(k)| 
to Fig.[5(i)]. 



D. Cooperative pearlite— martensite nucleation 

Below T = 0.5 martensite formation is followed 
by pearlite nucleation in a clearly sequential process. 
If xic = 1, at temperatures between T = 0.5 and 
T ^ 0.7 the nucleations of pearlite and martensite are 
simultaneous. In the case of X12 = xic = 1, at T = 0.6 
and t = 50 the system is mostly made of austenite, with 

At t 



a few pearlite and martensite nuclei. Fig. 12(a) 



100, the system is made of small pearlite and martensite 
grains. Fig. 12(b) By t = 200 austenite has disappeared. 
Fig. |12(c~ , and the microstructure is similar to what was 
obtained at lower temperature. Fig. ll(c)| but the path 
is noticeably different. At lower temperature martensite 
forms first then pearlite nucleates at a later time and at 
higher temperature pearlite forms first and martensite 
nucleates after the system is mostly pearlitic. 



VI. CONTINUOUS COOLING 
TRANSFORMATIONS 

So far, we have focused on isothermal transformations 
(TTT diagrams). However, these are not the most 
practical heat treatments. Instead of quenching the 
system and then holding the temperature constant, it is 
more natural to continuously decrease the temperature. 
If this procedure is repeated for various constant cooling 
rates, one obtains continuous cooling transformation 
(CCT) diagrams. The initial temperature is taken to be 
T = 1 in all simulations. As for TTT diagrams, we record 
and plot the times and temperatures corresponding 
to 10% martensite, 10% pearlite, and 10% austenite. 
Figure [13] shows the resulting CCT diagrams for several 
values of X12 and xic- 
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FIG. 13: Continuous cooling transformation diagrams for several values of Xi2 and xic- A: at least 10% austenite, M: at least 
10% martensite, and P: at least 10% pearlite. The thin lines correspond to the TTT diagrams of Fig. (2] 



Slow cooling rates, which correspond to a quasistatic 
situation, give rise to pearlite. On faster cooling, pearlite 
does not have time to form and martensite forms instead. 
One can notice that martensite forms at a temperature 
lower than the martensite start temperature of the TTT 
diagrams. This is because no martensite will form until T 
becomes smaller than 0.5. Then it takes a finite amount 
of time for martensite to form, during which temperature 
continues to decrease. This effect is more noticeable for 
very fast cooling and the discrepancy becomes smaller 
and smaller as the cooling rate decreases. 



to-pearlite lines tend asymptotically to their isothermal 
counterparts. 



The two dotted lines in Fig. |13(a)| correspond to the 
critical cooling rates required to obtain pure martensite 
and pure pearlite, respectively. At intermediate cooling 
rates, cooling is slow enough for pearlite to form but too 
fast to allow the transformation of austenite to pearlite 
to complete. At low temperature this retained austenite 
transforms to martensite, leading to the martensite- 
pearlite region seen in Fig. [I3(a)j One can notice that 
for long times, the austenite-to-martensite and austenite- 



For xi2 = 1 and xic = [Fig. |13(b)| , there is a 
martensite finish line, as was already observed in the case 
of isothermal transformations. One can also notice an 
austenite-pearlite-martensite mixture at intermediate 
cooling rates. No such mixture exists in the isothermal 
case. 



Figures 13(c) and 1 13(d)] show OCT diagrams for xic = 
1. Like the TTT diagrams shown in Figs. 2(c)| and |2(d)[ 
they exhibit several mixed microstructures. However, 
they are shifted to lower temperatures and longer times. 
For slow cooling, first pearlite forms, leading to an 
austenite-pearlite mixture. Then austenite disappears, 
leaving pure pearlite. Finally martensite forms, which 
gives rise to a pearlite-martensite mixture. On faster 
cooling, martensite and pearlite form cooperatively, as 
was already observed in the isothermal case. Soon after 
pearlite and martensite form, austenite disappears and 
the final microstructure is again pearlite plus martensite. 
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VII. CONCLUSION 

We presented a phase-field model which can be used 
to study alloys which can undergo displacive as well 
as diffusive transformations. It captures the important 
features of time-temperature-transformation diagrams, 
continuous cooling transformation diagrams, and mi- 
crostructures. It also sheds some light on the role of 
the interplay between the two types of transformations 
in stabilizing mixed microstructures. The existence of a 
martensite finish temperature (i.e. of retained austenite) 
is due to a hydrostatic strain associated with the 
martensitic transformation. When a strain is associated 
with pearlite formation, martensite and pearlite form 
cooperatively at intermediate temperatures, i.e. in the 
region of the TTT diagram where bainite is typically 



found in steel. The model also shows that in these mixed 
microstructures the habit planes are different from the 
pure martensite case and that small differences in volume 
changes can have noticeable effects on the early stages of 
martensite formation and on the resulting microstruc- 
tures. OCT diagrams show a shift of the transformations 
towards longer times, consistent with experiments. They 
can also exhibit mixed microstructures, which cannot 
exist in isothermal transformations. 
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